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Abstract. We show how to measure the order- two Renyi entropy of many-body 
states of spinful fermionic atoms in an optical lattice in equilibrium and non- 
equilibrium situations. The proposed scheme relies on the possibility to produce and 
couple two copies of the state under investigation, and to measure the occupation 
number in a site- and spin-resolved manner, e.g. with a quantum gas microscope. 
Such a protocol opens the possibility to measure entanglement and test a number of 
theoretical predictions, such as area laws and their corrections. As an illustration we 
discuss the interplay between thermal and entanglement entropy for a one dimensional 
Fermi-Hubbard model at finite temperature, and its possible measurement in an 
experiment using the present scheme. 
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1. Introduction 

Quantum degenerate gases, and ultracold atoms in optical lattices in particular, provide 
a unique framework to study quantum many-body physics [TJ El [3j H] . This refers first of 
all to the possibility of controlling many-body dynamics via external fields, thus allowing 
one to effectively engineer a wide class of interesting many-particle Hamiltonians, 
including those for strongly correlated systems [5j El [TJ [8] . Furthermore, a plethora of 
new measurement tools are available in atomic setups based on probing atoms with laser 
light, providing access to physical observables of many-body dynamics, in a way which 
is unparalleled in a condensed matter context. An outstanding example is the recent 
development of a "quantum gas microscope" for atoms in optical lattices, which allows 
single-atom detection and imaging with resolution of the lattice spacing in single-shot 
measurements [9], [10] - Given these unique and novel tools the challenge is now to identify 
new atomic measurement protocols that allow access to new many-particle observables 
of interest. Below we describe such a protocol, which allows the direct measurement of 
both thermal and entanglement Renyi entropies of (spinful) fermionic atoms in optical 
lattices for both equilibrium and non-equilibrium situations. The protocol consists of 
preparing two identical copies of the many-body systems in ID or 2D optical lattices, 
performing simple single particle operations that are readily implemented in optical 
lattices, followed by a read out with the quantum gas microscope. While in recent work 
[TT] we have discussed such a protocol for bosonic atoms, the fermionic case requires 
rather different arguments resulting in a different translation table to interpret the 
measurement results, although - quite remarkably - the basic procedure parallels the 
case of bosons. 

Direct measurement of thermal and entanglement entropy in atomic gases brings 
fundamental concepts, which so far have been discussed exclusively in a theoretical 
context, to the laboratory. Examples include the area law scaling that lies at the heart 
of the success of matrix product states methods [12] and logarithmic corrections in 
critical systems that allow us access to properties of the underlying conformal field 
theory [TBI [Hj. Also, the ability to explore higher dimensions is invaluable as the 
theoretical understanding of the corrections to the area law is less advanced and has 
sparked a lot of recent interest in the theory community [151 EE El [HI UHl EQl EH [22] . 
Other exciting applications include the possibility to detect topological order [231 El] or 
to monitor the dynamic generation of entanglement in quantum quenches [251 ESI E3 EE] 
- which are also relevant to questions concerning thermalisation in closed quantum 
systems [221 GHH US] - as already implemented in cold-atom experiments [32]- An 
important question in this context will be to ask which amount of the entropy is due to 
entanglement and what contribution is thermal entropy. This leads into the question as 
to whether it is possible to access the quantum entangled regime in current experimental 
setups, at least for small systems. 

To address the latter questions, we use quantum Monte Carlo (QMC) calculations 
to directly access the finite-temperature Renyi entropies for spinful fermions, described 
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by a Hubbard model in one dimension. By comparing the finite-T results to ground 
state DMRG calculations, we study the crossover between regimes dominated by 
quantum entanglement and regimes dominated by thermal entropy for realistic system 
parameters [33J . These simulation results allow us to give a critical assessment of when 
such a measurement could realistically be implemented by addressing the question of 
limitations of the measurement protocol. In particular, the exponential scaling of the 
number of single measurements with the entropy in the system sets the boundary 
in term of temperatures and system sizes, for what will be accessible under realistic 
circumstances. 

2. Measurement Protocol for Renyi Entropies of Fermionic Atoms 

We consider a many-body system represented by fermionic atoms in an optical lattice. 
The state of the system at a given time t is described by a density operator p. In the 
following discussion we leave the specific form of the state open. It can represent a pure 
state, p = \ip)(ip\, such as the ground state of a Hamiltonian, a thermal state, or any 
other (mixed) non-equilibrium state of the fermionic atoms. Our goal is to develop a 
protocol to measure Renyi entropies both for the total system and for subsystems. For 
the total system the Renyi entropy of order two is defined by S^p) = — logTr{p 2 }, 
and thus is given in terms of the purity of the density operator P2G ) — Tr{p 2 }. For a 
subsystem 7Z we define a reduced density operator p-ji = Tr^.{p}, and a corresponding 
Renyi entropy as S^pn) — — l°g Tr^jp^}. The knowledge of both S 2 (p) and ^(p^), 
allows one to quantify the entanglement of the subsystem 7Z with the rest [SHI35J. While, 
as emphasized above, the following discussion is valid for any quantum state, specific 
scenarios of experimental interest include monitoring state purity and entanglement 
entropies in quench dynamics as a function of time, or thermal vs. entanglement entropy 
in thermodynamic equilibrium situations. 

The Renyi entropy and purity are nonlinear functionals of the quantum state and 
thus not directly observable. However, the purity can be directly obtained by measuring 
several copies of the same state [35] . It can be expressed as the expectation value of 
the swap operator V2 on a system that is prepared in two identical copies in the same 
quantum state p, that is, Tr{p 2 } =Tr{V2p <E> p} = (V2), where the swap operator is 
defined as V^i) ® ^2) — IV^) ® Similarly, the purity of the reduced density 

operator of a subsystem TZ is given by the expectation value of the operator that 
swaps the quantum states just in the part TZ, that is Tr{p^} = TrjV^p ® p} = (V^). 
Measurement of Renyi entropies S^p) and ^(p^) of fermionic atoms in optical lattices 
thus reduces to (i) the ability to prepare two identical copies of the atomic system, and 
(ii) the development of protocol to measure (V2) or (V^) by simple operations and read 
out in an optical lattice on the two copies, in a way that can be readily implemented with 
high fidelity in present experiments. While the first aspect is primarily an experimental 
question, we will focus in the following on the protocol to realize the measurement of 
the swap operation (ii) and interpretation of measurement results of a quantum gas 
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microscope to determine (V2) or (V^ 1 ). We comment on the assumption of identical 
copies (i) at the end of this section. 

To be more specific, we consider spin- 1/2 fermions in an optical lattice, in a setup 
illustrated for two 2D lattices representing the two copies in Figure [IJ We define a basis 
of Fock states for the two copies of the system, 

|n,m}^nK^n( & U mj ' CT | vac )' C 1 ) 

where a\ a and b\ a denote the creation operations for fermionic atoms on lattice site i 
and spin a = ±1/2 in the first and second system (copy), respectively, and with the 

n = {7ii,f, 711,4., ri2,t> • • •} an d m = { m i,t> m i,4-> m 2,t; • • •} are the occupation numbers of 
the two systems. The swap operator V2 acts on these states as V^n, m) = |m, n), which 
interchanges the configurations n and m. For bosons this swap operation amounts 
to the interchange a\ a <H- b\ a . Ref. [TT] has shown that this is readily implemented 
in an optical lattice by turning on tunnelling between corresponding lattice sites of 
the two copies, and reading out lattice occupation of the first copy (modulo 2) with 
the quantum gas microscope. However, for fermions applying V2 is not equivalent to 
exchanging a ia -H- b ia , since there is an ordering problem, i.e. one must keep track of 
the fermionic signs. In the following we will present this protocol for fermions. We will 
present our results first in the form of an experimental recipe in Section 2.1 and give 
the formal proof in Section 2.2 and Appendix A| In Section 2.3 will analyze limitations 
and scaling of errors in these measurements. 



2.1. Experimental Protocol to Measure Tr{p^} = Tr\V^p £g> p} 

The expectation value of is obtained by averaging over a series of single 
measurements, where each single measurement proceeds in three steps as illustrated 
in Figure [TJ 

(i) Initially two identical copies of the same low entropy (non-) equilibrium state 
are prepared [37]. With optical lattices this can be performed in two parallel ID tubes 
or two neighbouring 2D planes that are completely decoupled, (c.f. Figure [l|a). 

(ii) At a given time t the lattice depth within each copy is suddenly ramped 
up to freeze the atomic configuration by suppressing tunnelling (atomic limit). 
Simultaneously, all interactions (e.g. between different spins) are turned off, e.g. by 
using magnetic or optical Feshbach resonances [3"S] I3"U] . Alternatively the protocol can 
be executed on a timescale where the effects of interactions are negligible. Note that 
during this step the entanglement and the overall entropy does not change. 

We then lower the barrier between the two copies to allow tunnelling between each 
site and its copy with an amplitude J a b for a fixed time T a b = 7r/(4J a f,) (c.f. Figure [Tjb) . 
This can be accomplished with the use of an optical superlattice [ID]. With this 
operation we implement the beam splitter U2 that maps 

U 2 ■ a ii<T ->■ A=(a it(T + b ita ); b i>a ->■ -j=(b i)ff - a ii(J ), (2) 
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Figure 1. (color online) The three steps in the measurement protocol of second Renyi 
entropy for spinful fermions in a 2D optical lattice, (a) First, two copies of the many- 
body state are produced, (b) Then the barrier between the copies lowered such that 
atoms can tunnel from each mode to its copy to realize a beam splitter operation 
between the copies (see Equation ([2])). (c) Finally, the fermion number is measured 
site- and spin-resolved in both copies with a quantum microscope. Here we show a 
typical outcome of a single measurement run, where on each site on finds either no 
atom, an atom with spin up, an atom with spin down or both spins on one site. In 
this example, according to Table[T]the measurement outcome for the swap operator on 
the whole system is +1, while the result for the reduced set of modes that are enclosed 
by the dashed line is —1. Simultaneously one also obtains measurement results for all 
other subsets. For example the outcome for the swap operator on the spin up modes 
is -1. 



where aj jCT denotes the annihilation operator at site i with spin a and bi yCr denotes the 
annihilation operator in the corresponding mode in the second copy of the systemTJ] 

(iii) Finally we measure site- and spin-resolved the occupation numbers using a 
quantum gas microscope [9j [10]. We denote the measurement results in copy one by 

In the case of fermions this number is either zero or one 



n\ a and in copy two by n\ 



(2) 

2,(7 * 



X To strictly realise (|2j) one needs additional phase shifts. They can be obtained by shifting the energies 
of the two copies relative to each other in an additional step, e.g. using a superlattice. However, due 



to the particle number superselection rule (see Appendix A) these phase shifts do not affect the final 



measurement result, and are not necessary for the present protocol. 



A Measurement Protocol for Entropies of Fermionic Atoms in Optical Lattices 



6 



due to the Pauli principle, but also for bosons it suffices to determine the parity of the 
occupation number. 

The difference between bosons and fermions consists of how the measured 
occupation numbers relate to the measured value of the operators V 2 n . 

Bosons. If one finds an even (odd) number of bosons on the modes belonging to TZ of 
copy one, then the measurement outcome for V 2 n is plus (minus) one [11] . 

Fermions. For fermions the measurement outcome of the operator V 2 n depends on the 
total number of fermions in both copies of the modes belonging to TZ, which we denote 
by N^ t = a)en n i)r + ra i>' an d on f ne number of atoms on copy one of the modes 
belonging to TZ, denoted by Nf- = a ^ n n^. The corresponding measurement result 
for V 2 TL can be read off from table fil 



iv tot 




N? 


result for V 2 n 


even 


even 


even 


1 


even 


even 


odd 


-1 


even 


odd 


even 


-1 


even 


odd 


odd 


I 


odd 










Table 1. Rules determining the measurement outcome for fermions. For a derivation 
see Section [2T2l 

To determine the expectation value of V 2 n , and thus the purity of the corresponding 
reduced density operator, one has to repeat the whole measurement procedure and 
average over the outcomes. Since all the swap operators for different subsets TZ commute 
([V^, = 0), it is in principle possible to measure all of them at once. In fact in 
each single run one obtains a measurement result for all possible subsets. 

2.2. Justification of the measurement protocol 

Here we show that the above protocol implements a measurement of the swap operator 
and thus the Renyi entropy. We give the proof for V2 and at the end comment on V 2 n . 

As we noted above, the swap operator V 2 acts on the Fock states according to 
I4|n, m) = |m, n). In the bosonic case it is simply given by the operator which 
interchanges a\ <H- b\, since the order of the creation operators does not matter for 
bosons. However for fermions applying V 2 is not equivalent to exchanging a\ <-> b\, 
since such an operation performs the mapping |n, m) — > (— 1)^ 4 j |m, n). 

From the Fock-basis one can easily construct the eigenbasis of V%. The eigenspace 
with eigenvalue +1 is spanned by vectors \tp^ m ) which are of the form |^ n ) = |n, n), 
and iV^m) = r^(l n ' m ) + l m ) n )) f° r n 7^ m. The eigenspace with eigenvalue —I 



A Measurement Protocol for Entropies of Fermionic Atoms in Optical Lattices 



7 



is spanned by the vectors \ipn,m) = ^(|n, m) — |m, n)). Under the beam-splitter 
operation Equation §2§ the eigenstates of V 2 transform into a superposition of states 
in the occupation number basis according to 

W£ J = E C £n,ml k ' ^ + HI - k). (3) 
k 

The coefficients c^. n m depend on whether we are considering bosons or fermions. We 
discuss the two cases separately. 

Bosons. The coefficients in Equation ^ are given by 

^n,m=(l±(-l)^*0^m, (4) 

where G?k ; n,m is a numerical factor which is irrelevant for the following discussion. From 
equations Equation ^ and Equation Q we see that the beam splitter operation (BS) 
transforms the symmetric states |^n m ), that is the eigenspace with eigenvalue +1 into 
the space with an even number of atoms in the first copy, since all coefficients for states 
with an odd number of atoms in copy one, Yli vanish after the application of the BS. 
Similarly the eigenspace with eigenvalue —1 is transformed into the space with an odd 
number of atoms in the first copy. 
Fermions. We find for fermions 

Ck;n, m = (1 ± (-l)S* e k;n , m , (5) 

where ek ;n ,m; as for bosons, is a numerical factor which is irrelevant for our purpose. 
The main difference to the bosonic case is the dependence on the parity of the term 
^ijUiTrij in the exponent. If this term is even, the situation is the same as for bosons, 
and the symmetric states are mapped onto states with an even number of particles in 
copy one, while antisymmetric states are mapped onto states where this number is odd. 
However if ^ . nirrij is odd, the situation is reversed. Unfortunately there is no way of 
determining the parity of Ylij n i m j after the beam splitter operation has been applied. 
One only has access to the values of ki + U and rij + m^, that is the number of particles 
in the modes of copy one after the beam splitter, and the total number of atoms in both 
copies of each mode (which is the same before and after the beam splitter operation). 
However, we will show in the following that this sufficient information to proceed. The 
basic idea is the following. 

(i) First consider only the eigenstates |^ m ) with Ylii n i = 12i m i = ^ that is 
with an equal number of atoms in the two systems (before the beam splitter operation 
is applied). For these eigenstates the parity of Ylij n % m j — N 2 can be determined from 
the knowledge of the total number of atoms N to t = 2iV, which can be accessed from the 
occupation number measurement in the final step of our protocol. If N is even, then 
N 2 is even as well, and as in the bosonic case, the (anti-) symmetric states are mapped 
onto the space with (odd) even number of particles in copy one. On the other hand, 
if N is odd, then also N 2 is odd, and the situation is reversed, the (anti-) symmetric 
states are mapped onto the space with (even) odd number of particles in copy one. This 
is reflected in the rules presented in table [TJ 
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(ii) Second, consider only the eigenstates |"0^ m ) w hh ^2i n i ^ X^ m i- Then the 
total number of fermions N tot is not enough to determine the parity of Ylij n i m j- 
Thus, when assigning a measurement outcome according to table [TJ there are some 
pairs (n, m) where one incorrectly assigns a measurement result of +1 to the state 
l^n,m) an d a value of —1 to the corresponding state |V>n,m)- However, for product 
states constrained by a particle number superselection rule [JT] it is readily shown (see 



Appendix A) that the probability of finding the system in the symmetric state \ipnm) 



is the same as the probability of finding it in the corresponding antisymmetric state 
iV'nm)) ^ Ylii n i 7^ Ylii m i- Therefore these instances average to zero and the error is 
irrelevant for the average value. Furthermore, if the total number of atoms iV tot is odd, 
then one can be sure that £V rii ^ ^ rrii, and assign a measurement value of zero right 
away, as suggested in table [1} A more formal and detailed proof of these points can be 
found in Appendix A In contrast to the bosonic case [llj, the above arguments can not 
be generalised in a straightforward way to higher order Renyi entropies for fermions. 

If one is interested in the purity of the reduced density operator of a subset of 
modes TZ, that is in the measurement of V^ 1 , the discussion is completely analogous. 
Note that, since the the beam splitter is unitary and the beam splitting operations for 
different modes commute, it is irrelevant whether or not one performs the beam splitter 
operation also on modes not belonging to TZ. Thus in practice one can perform always 
the full beam splitter U2 on each pair of modes, and determine the number of atoms and 
the corresponding measurement result for V^ 1 in each subdivision TZ simultaneously. 



2. 3. Limitations and Effect of Errors 

The main limitation of the proposed measurement scheme is that the number of single 
measurements necessary to determine the entropy with a certain statistical accuracy can 
become large. This number can be estimated as follows. Each single measurement gives, 
according to table [l] either ±1 or 0. Their mean value determines (V^ 1 ) = Tr{p^}. The 
variance of the measurement outcomes can be expressed as (AV^) 2 = (P evcn ) — (V] 1 ) 2 , 
where P ev en is the projector on the sector with even total number of particles in the two 
copies. It is easily shown for product states of two copies, both constrained by a particle 
number superselection rule, that 1/2 < (P CVCQ ) < 1. To determine the expectation value 
of and thus the purity with a certain statistical accuracy ay, one needs a number 
of measurements j^y that is given by j^ydy = (P eve n) / (V^ 1 ) 2 — 1. For highly mixed 
states this number diverges as j^yo v ~ 1/Tr{p 2 }. This results in an exponential scaling 
of the required number of measurements with the Renyi entropy. To determine the 
entropy with a relative statistical uncertainty of as one needs #5 measurements with 
#sO"i = ^2 ((-Pevcn )e 252 — 1) ~ e 25 2 -2iogS 2 _ ^y e note fast all known schemes to measure 
entropy based on multiple copies [36l [TT| W2\ suffer from this limitation. 

Since the (parity of the) number of atoms on each site enters crucially in the 
measurement result, errors in the measurement of this numbers are a major error source 
in an experiment. Their effect is most easily outlined in the measurement scheme 
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for bosons, but the discussion can be carried out analogously for fermions. In an 
ideal experimental implementation the system state determines the probability p± of 
finding an eigenvalue plus or minus one, such that (V2) = p+ — p~. Suppose that 
with a probability e the quantum gas microscope incorrectly measures an even (odd) 
number of particles on a certain site, when there is actually an odd (even) number 
of atoms. Assuming that such errors occur in an uncorrelated fashion and with the 
same probability on the M sites that are measured, the experimentally determined 
expectation value (V^e is reduced compared to the actual one by (V^e = (V2HI — 2e) A/ . 
Thus the measured entropy S^eO ) is just the sum of the actual entropy of the 
system S 2 (p) and a contribution from the quantum gas microscope S'niicroscope = 
— M log(l — 2e), where the contribution of the quantum gas microscope is extensive in 
the size of the measured (sub) system. Thus, the measured purity (entropy) is always 
smaller (larger) than the actual one. 

A general assumption [13] underlying the protocol is the preparation of two perfect 
copies p ® p. Even though our protocol does not allow to check this assumption, it 
can strengthen it a posteriori, e.g. if the two copies are close to a pure state. Then 
the measured expectation value is one if the state in both copies is the same. Also, 
one can relax this assumption and still obtain useful information from the described 
measurement protocol as outlined in Appendix B On the other hand, one might be 
interested in the overlap of two states that are different from the outset. By preparing 
the two states in the form p\®p2, this protocol gives access to Tr{p!p 2 } = Tr{V2Pi<S>P2}- 



3. Renyi Entropies of a Fermionic Hubbard Chain 

The protocol to measure Renyi entropies for fermionic systems described in the previous 
section provides a novel tool to access in experiments fundamental properties of many- 
body systems related to entanglement and thermal entropy. One of the important results 
quantum information theory has brought to the field of quantum many-body systems 
is the finding that ground states of local Hamiltonians typically exhibit an entropic 
area law S ~ adA, where OA denotes the perimeter of the boundary delimiting the 
two complementary subregions A and B |44j . In the particular case of one-dimensional 
critical systems admitting a conformal field theory (CFT) description, the area-law 
picks up an (additive) logarithmic correction, whose prefactor solely depends on the 
central charge c of the CFT. These are very important theoretical results underlying 
the success of matrix- and tensor network based numerical and conceptional methods, 
and have furthermore deep connections to quantum field theory, string theory and 
black hole physics. It is therefore highly desirable to test these theoretical predictions 
in actual experiments, enabled through the experimental protocols introduced in 
previous [361 M H2] and the present work. 

In the present section we illustrate this for the example of the Fermi-Hubbard model 
in one dimension, as it is the simplest model of interacting fermions that can be realised 
with cold fermionic atoms in an optical lattice. In terms of creation (a ia ), annihilation 
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(oi,(r) and counting operators (n^, 
lattice with L sites is given by 



a ia ai ja ) the Fermi- Hubbard Hamiltonian for a ID 



t 



L 



"=-«'EE(' 



a| )0 .ai + i )CT + h.c.j + C/ y^n^nj,;, 



(6) 



a=t,4- »=1 



i=l 



where tp denotes the hopping amplitude between neighbouring sites and U the onsite 
interaction energy. The model is exactly solvable via the Bethe ansatz [451 [46] and the 
phase diagram exhibits a metallic two-channel Luttinger liquid ground state at generic 
fillings for all interactions U > [47, 46J. At half filling and U > the charge degrees 
of freedom are gapped and a Mott insulator appears. 

Here, we use a generalised directed loop algorithm within the stochastic series 
expansion (SSE) framework [injHniED] to access the thermal Renyi entropies following a 
measurement scheme based on the dynamic update of the world line topologies presented 
in Ref. [21] . Large blocks are built up consecutively using the increment trick [51] that 
allows for an efficient update in replica space. 

In the following we quantify the ground-state entanglement as well as the 
thermal Renyi entropy that one would obtain with the proposed protocol for realistic 
experimental system sizes. We show to what extent quantum entanglement can be 
accessed through a measurement of S2 for system sizes and temperatures available in 
an experiment. A systematic study of the crossover between entanglement and thermal 
entropy is presented in [33J. 

3.1. Zero Temperature 

Let us start the by looking at the Hubbard chain at zero temperature. Figure [2] shows 
the n = 2 Renyi profiles for bipartitions A U B with block sizes I a obtained using 
DMRG for chains with L = 48 sites, both for half (n-f — — 1/2) and quarter filling 
(n-f- = — 1/4). A prominent feature of the the Renyi profiles for finite system sizes 
is that they exhibit characteristic oscillation associated with the Fermi- momentum kp, 
giving rise to two (four) branches at half (quarter) filling. Further, one clearly identifies 
the envelope carrying the logarithmic corrections to the area law. 

These features are well understood via the underlying conformal field theory 
[52] [53] [54J [55], from which the Renyi profiles are obtained to be 



where l' A = L/tc sin(7rZ a/ L) is the chord distance. The first term is the leading 
contribution and describes the logarithmic increase of the entropy with block size. It 
is directly proportional to the central charge c and distinguishes between the metal 
(c = 2) and the Mott insulating state (c = 1) having two respectively one gapless 
channels. By measuring the Renyi profiles with the previously introduced protocol it is 
in principle possible to determine this pre-factor and to extract the central charge c in 
an experiment. This is conveniently done by plotting S2 as a function of log(2/^), as 
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half filling quarter filling 




Figure 2. The upper left (right) panel shows the second Renyi entropy 52 for the half- 
filled (quarter-filled) Hubbard model as a function of the chord distance l' A for L = 48 
and U/tp = 8 obtained using DMRG. The lower panel displays the same entropy data 
as a function of \ogl' A , where l' A is the chord length of block A. Here, we show DMRG 
data for systems up to L = 96. The slope of S 2 is given by c/8 (see Equation Q) 
and the blue and green lines are guides to the eye corresponding to c = 1 and c = 2 
respectively. 



shown in Figure [2j This way, 52 approaches a straight line with slope c/12(l + 1/rt) for 
large block and system sizes. One can see from Figure [2] that the data for the half filled 
Hubbard Model at U/tp = 8 is consistent with a central charge of c = 1, as the two 
branches of the Renyi entropy approach the corresponding asymptotic line from above 
and below whereas the quarter filled chain exhibits c = 2. 

3.2. Finite Temperature and Trap 

A measurement of the Renyi Profile as proposed in Section [2] will necessarily be 
performed at finite temperature. Then, the Renyi entropy will not only pick up quantum 
entanglement but also thermal contributions to the entropy. Since the thermal entropy 
is an extensive quantity, one expects a crossover from the area-law to a volume law. This 
is shown in Figure |3j where we plot finite-temperature QMC results for U ftp = 6. To 
make a closer connection to the experimental situation, we include a harmonic trapping 
potential, i^trap = Vi n u with V(i) = (i — L/2) 2 /20 — 6. Filling the trap with 18 



A Measurement Protocol for Entropies of Fermionic Atoms in Optical Lattices 



12 








8 



16 



24 



32 



Figure 3. upper panel: Renyi profile for a Hubbard chain in the presence of a trapping 
potential with V(i) = [i — L/2) 2 /20 — 6 and U/tp = 6. The green crosses correspond 
the ground-state of the system obtained from DMRG with 18 particles whereas the 
finite-temperature results are obtained from QMC using grand-canonical simulations 
with an average particle number coinciding to the T = DMRG calculations. A 
density profile for different temperatures is shown in the inset, lower panel: Mutual 
information for two disjoint blocks aligned with the edges of the system for different 
temperatures as depicted in the illustration. 



particles, at zero temperature a Mott insulator extends over about eight lattices sites 
and the density drops to zero in the wings quite rapidly (see the density profile in the 
inset of Figure [3]). Correspondingly, the T = Renyi profile shows the characteristic 
2kp oscillations in the centre of the trap that indicate the proliferation of the Mott 
phase. The metallic wings that appear due to the presence of the harmonic trap, are 
also clearly visible in the Renyi profile in Figure [3} 

When the temperature is increased above the finite size gap, the entropy S2 picks 
up a linear contribution (see Figure [3J. Thus, the 5*2 profile is no longer symmetric with 
respect to the centre of the system and the purity of the whole system decreases as S^L) 
increases. In addition the amplitude of the parity induced oscillations is suppressed as 
the temperature is increased. One can see that for temperature T/tp > 0.1, the entropy 
is already dominated by a linear increase of the thermal entropy and 5*2 looses its parity 
effects. 

The ability to measure Renyi profiles also allows access to other quantities that 
provide further insight into the entanglement and correlation properties of the system. 
For instance, the mutual information between two (possibly disjoint) blocks A and B 
and is given in terms of the Renyi entropies as 



This measure is particularly useful when dealing with mixed states as it does not pick 




(8) 
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Figure 4. (left) Renyi entropy of the whole chain, S%(L) at quarter filling and 
U/tp — 4, for different system sizes L at different temperatures T. The contour lines 
denote constant temperature T/tp = 0.1, 0.15, ... 0.5. (right) Number of measurements, 
#5, required for a measurement of Sa with a statistical relative error as, using our 
protocol. 



up an extensive contribution from the thermal entropy - it obeys an area law [561 EH HI] 
even at finite-T - but is sensitive towards correlations between the two subblocks [58J. 
For a detailed discussion of the mutual information in the Hubbard chain we refer to 
Ref. [33]. 

3.3. Limitations 



As pointed out in Section |2.3[ the ability to determine these entropy profiles using the 
protocol proposed in this work is fundamentally limited by an exponential growth with 
entropy of the number of single measurements required to obtain a certain statistical 
accuracy as (see Figure To give numbers, for example 10 5 single measurements 
are needed to determine Renyi entropies up to 52 ~ 5 with a relative statistical error 
of as ~ 0.1. For one dimensional systems, at temperatures below the finite size gap 
this is not a sever restriction as the entropy obeys a area law with at most logarithmic 
corrections. For example in the previous subsection we showed that in the ID Hubbard 
model, at temperatures below the finite size gap the entropies are typically of the order 
one. This values can be resolved with just ~ 10/a| single measurements. Above 
the finite size gap however the number of measurements sets the limit in terms of 
temperatures and system sizes, as the entropy becomes extensive in the system size and 
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increases with temperature. 

To quantify this, Figure [4] shows the Renyi entropy S 2 (L) of a full (homogenous) 
Hubbard chain at quarter filling and U/tp = 4 as a function of system size L (open 
boundary conditions) and temperature T, as well as the number of measurements, #5, 
required to resolve a certain value of S 2 with a statistical relative error, as, using the 
above proposed protocol. This clearly shows that experimentally relevant situations, 
e.g. T /t F ~ 0.2, L ~ 30 can be explored with a moderate number #5«t| < 10 3 of single 
measurements. 

4. Discussion 

To summarise, we have presented a method to determine the order two Renyi Entropy 
for bosons and fermions in an optical lattice. The scheme is based on the possibility of 
preparing two identical copies of a quantum state in an optical lattice, coupling the two 
copies via a superlattice, and site-(and spin-) resolved measurement of the occupation 
number (modulo two). The combination of these tools, which are available in current 
experiments, allows to directly determine the entropy of a many-body quantum state, 
and may be seen as a thermometer for states near the absolute ground state. On 
the other hand, for pure states this opens the possibility to study entanglement as 
quantified by the entropy of a subsystem. Possible applications include test area laws 
for the scaling of entanglement entropy in ground states, and monitoring entanglement 
in time-dependent systems. Using QMC-techniques we analyzed the possibility to test 
area laws and their corrections in a finite size system at finite temperatures typically 
present in current experiments. 
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Appendix A. Details on the Proof of the Measurement Protocol 

Here we elaborate more formally on some points of the proof presented in Section |2.2| 
for fermions. First, note that a pure state (in one copy) with N fermions (constrained 
by the particle number superselection rule) always has the form: 

ri 
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where N is the total number of atoms. A general state is a mixture of such states (with 
possibly different total number of fermions): 

p = H E <°2>>H- ( A - 2 ) 

N n, m 

E* «> = N 

J2i m i = 

Thus, a general product state in the two systems has the form 

P®P=E E E piM£hfi>Rm|. ( A -3) 

JV,iV' n, m n, m 

Ei «i = w Ei «» = $ 

Ei m i = W l»i = Jv 

They have the crucial property that for Yli n i 7^ Si m « one has 

(^n.mlP ® PlV'Jm) = (^n,mlP ® P|^n,m)> ( A - 4 ) 

as one can easily show using the definition if the eigenstates (see Section 2.2[ ) and 



Equation (A.3). This property is used to allow for errors in the assignments of the 
measured eigenvalue, that do not alter the average value. An alternative way of thinking 
about this is the following. We want to determine the purity of p via an the expectation 
value Tr{p 2 } = Tr{V2p <S> p}. However V 2 is not the only operator whose expectation 
value in the state p <g> p is equal to Tr{p 2 }. Below we will introduce a whole family of 
operators V 2 (depending on some function /), which share this property. We define 

P±(n, m) = -(|n, m) ± |m, n))((n, m| ± (m, n|), (A. 5) 

and further Q(n,m) = |(P + (n, m) — P_(n, m)), such that we can write V 2 = 
^2 nm Q(n,m). The factor 1/2 in the definition of Q is to correct for double-counting. 
We can divide the eigenstates of V 2 in two classes, those with £\ n i = J2i m i an d the 
others, leading to the representation 

v 2 = Q(^ m ) + E Q(^ m )- ( A - 6 ) 
11 in n in 

Ei n i = E« m i Ei n i i 1 Ei m i 



Note that due to property Equation (A.4), the expectation value of each term in the 



second sum is identically zero for states of the form Equation (A.3). Therefore, the 



expectation value of the Operator V 2 is the same as the one of Vj 17 defined as 

V 2 if) = J2 Q(^™)+ E /(n,m)Q(n,m), (A.7) 
11 in 11 in 

E; »> = Ei m t Ei i 1 Ei m i 

where /(n, m) is an arbitrary (real- valued) function. In particular we can choose 

s(n + m) = 1 

f(n,m)={ 1 s(n) = s(m) = s((n + m)/2) , (A.8) 
-1 s(n) = s(m) 7^ s((n + m)/2) 

where we use the notation s(x) = SjXj mod 2. The eigenstates of V 2 are the same 
as those of V 2 , but they have different eigenvalues. The advantage of choosing / in the 
above form is that the eigenvalue of V 2 can be determined form the total number of 
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atoms and the number of atoms in copy one after the beam splitter, as presented in 

table [Tj This can easily be checked by explicitly looking at the transformation of all the 

( f) 

different classes of eigenstates of V 2 under the beamsplitter. 



Appendix B. Mixtures of Copies 

As pointed out in [36] the expectation value of the operator Vi in the state p to t = p £S> p 
is given by Tr{p 2 }. Although this assumption is standard it is difficult validate [43J. 
For example fluctuations in the intensities of the laser beam that generates the lattice 
[59] affect the two copies in the same way and lead to correlated errors in the state. 
Such identical, but correlated errors in general lead to a state which is a mixture of 
different products of the same state, that is p to t = J2iPiPi ® Pi with YliiVi — 1- I* 1 
the following we show that even in this situation the measurement of the operator 
provides useful information about the entropy in the sense that it provides bounds 
on the purity of the reduced system of one of the two "copies", that is of the state 
p = Tr copy2 {p tot } = YjiViPi- To this end we show that 

^Tr{U 2 p tot } < Tr{p 2 } < Tr{V 2 p tot }. (B.l) 
We first show the upper bound. Consider the function 

F({r}) = Tr{p 2 } - Tr{V 2Ptot } = J2p iPj Tr{ PiPj } - ^ ft Tr{p 2 }, (B.2) 

i,j i 

defined in the polytope spanned by the pi with YliPi = 1. Note that this function is 
zero at the corners of this polytope given by pi = 5y. Further, the hessian matrix of F 
is constant and given by iiy = d ^ = 2Tr{piPj}. This is a Gramian matrix and thus 
positive semidefinite. Thus the function F is convex everywhere and its value on the 
(convex) polytope defined by ^2 i p% = 1 is smaller than its value at the corners. This 
proves the upper bound. To prove the lower bound we calculate the (unique) global 
minimum of F. It is assumed at the position that satisfies dF/dpi = ^ . 2pjTr{piPj} — 
Tr{p 2 } = 0. That is, at the minimum we find pf£t{ P i P j} = ^Tr{p 2 }. Using 



this in Equation (B.2) we find that F > — | X^P«Tr{p 2 } = — |Tr{V2Ptot}, and thus 
|Tr{T^p tot } < Tr{p 2 }. These inequalities trivially hold also for purities of subsystems 
of one of the two copies (with the corresponding swap-operators). 



The inequalities Equation ( |B.l ) give bounds in terms of the expectation value of 



V2 on arbitrary mixtures of copies. In the typical experiment situation one expects to 
be very close to a simple product of two states with a small admixture of other copies, 
such that p = 1 — e with e = J2i^oPi ^ 1- Because the function F({pi}) is convex we 
can bound 

F({ Pi }) > F(Ro}) + - ko)^r^- (B.3) 

= - ^ Pi Tr{(p, - po) 2 } > -2 5>, = -2e, (B.4) 
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and thus one has Tr{\/ 2 p tot } - 2e < Tr{p 2 } < Ti{V 2 p tot }. 
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